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We discuss the procedure of extracting the photoproduction cross section for neu¬ 
tral pseudoscalar mesons off neutrons from deuteron data. The main statement is 
that the final-state interaction (FSI) corrections for the proton and neutron target 
are in general not equal, but for vr^ production there are special cases were they have 
to be identical and there are large regions in the parameter space of incident photon 
energy and pion polar angle, 9*, where they happen to be quite similar. The correc¬ 
tions for both target nucleons are practically identical for production in the energy 
range of the A(1232)3/2'’' resonance due to the specific isospin structure of this exci¬ 
tation. Also above the A-isobar range large differences between proton and neutron 
correction factors are only predicted for extreme forward angles {9* < 20°), but the 
results are similar for larger angles. The case of rj photoproduction is also shortly 
considered. Numerical results for the 'yp —>• ir^p and yn correction factors 

are discussed. Also the model description for the available data on the differential 
'yd —)• ir^pn cross sections are given. 


1. Introduction 

We analyze the procedure of extracting 7r°/?7-photoproduction cross sections off neu¬ 
trons and off protons from deuteron data. In both cases, this procedure takes into account 
the fact that the elementary yiV reaction can occur on the proton as well as on the neu¬ 
tron of the deuteron, and that there are effects from the hnal-state interaction (FSI). Thus, 
even in quasi-free kinematics (close to the kinematics of the process on the free nucleon), 
there can be corrections due to this nuclear effects. In the case of photoproduction off 
the proton, one can compare the cross sections for free protons to that extracted from 
the reaction off quasi-free protons bound in the deuteron. Thus, in the proton case, there 
is a possibility to verify the calculation of the corrections. This verihcation procedure is 
impossible for the reaction off the neutron, since free neutron targets do not exist. The 
strategy for the extraction of free-neutron cross-section data from quasi-free reactions off 
the deuteron is then to test the FSI modelling for the proton case and apply the same 
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Figure 1: The impulse-approximation (Md, Ma 2 ), A^A^-FSI {Mb) and ttA^-FSI (Md, Mc 2 ) 
diagrams for the reaction 7 p —)■ n^p. See text for dehnition of “fast” {Nf) and “slow” {Ng) 
hnal-state nucleons. 

model to the neutron. The aim of the present note is to employ a specihc approach to the 
FSI corrections for the deuteron data to obtain and compare the correction factors for the 
elementary photoproduction processes on the proton and neutron targets in wide regions 
of initial photon energy and of outgoing-pion angles. 

2. Reaction Amplitude 

Let us write down the ■yd —)• n^pn amplitude (the case of p production will be discussed 
later on) as 

Tf'yd = Mai + Ma2 + Mb Mcl Mc2-i (1) 

where the terms in the r.h.s. are represented in Fig. [H The amplitude M^d contains 
the impulse-approximation (lA) [1] terms Mai and Ma2 and FSI terms, where Mb is the 
A/V-rescattering, while Md and M^ are the vrA-rescattering amplitudes. Hereafter, we 
consider the kinematics with fast and slow hnal state nucleons Aj(qr^) and As(qr^), where 
Qjg are their momenta in the laboratory system, and \qj\ 3> Ig^l. First, let us discuss the 
elementary process yN -y ttN . In the isospin basis, the amplitude has the form [2] 

A{yN -y TlaN) = A^6a3 + ^vl 2 [A; Ts] + , (2) 

where is the Pauli matrix; Ay and Ayi are two isovector amplitudes and Ag is the 
isoscalar one. For the proton and neutron targets, we have 

A{yp 7r°p) = Ay + A^, A{yn -y 7r°n) = Ay - A^, (3) 

i.e., in the general case, when A^ ^ 0, the yp and yn amplitudes are not equal. Now, we 
shall consider two different kinematics. Case 1 corresponds to a fast proton and a slow 
neutron in the hnal state, i.e., Nf =p and Ng = n. We dehne case 2 as the result of the 
isospin replacement p n for case 1, i.e., Nf = ri and Ng=p, with identical momenta and 
spins of the fast and slow nucleons. Then, the lA amplitudes Mai and Ma2 read 

Mai = + klW) (fdiqg)], Ma 2 = [(A(,2) - A(2)) (fdiqf)] (case 1); 

(4) 

Mai = - A(b) (Pd{qg)], Ma 2 = [(A(,2) + A(2)) Lpd{qf)] (case 2). 
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where y:>{q) is the deuteron wave function (DWF). Hereafter in this paper, we consider only 
the isospin structure of the amplitudes and use the short notation [•••], which means that 
all particle momenta and spins are properly taken into account. More detailed expressions 
with all the variables written out can be found, e.g., in Refs. |3HS] (see also references 
therein). The superscripts “(1)” and “(2)” of the terms and mean that they are 
taken at different values of spin and momentum variables due to the replacement Nf Ng, 
i.e., A^J ^ A^^g. Both amplitudes Mai and Ma 2 change from case 1 to case 2 because the 
relative sign of the term Ag changes. Since |qr^| \Qf\y we have (p{qg) 3> ^{qf) and 
|Mai| S> |Ma 2 |. Thus, the leading diagram is Mai, i-e., the reaction yd —)■ ir^pn proceeds 
mainly through the subprocess yp —?• vr^p or yn —)■ vr^n in case 1 or 2, respectively (the 
detectable recoil nucleon is the participant of the reaction). 

Now consider the NN-FSI diagram M^, in Fig. [H First, let us introduce the NN- 
scattering amplitude in the form 

M{NiN2 ^ NsN,) = ^Mo{xtxl){X2'^Xi) + ^Mi(x3+TxD(x^2+rxi). (5) 

Here: Xi _4 are the nucleon isospinors; = T 2 X*, t = (h, 'r 2 , t's), where ri _3 are the isospin 
matrices; Mq and Mi are the isoscalar and isovector NN amplitudes. Making use of Eqs. (E]) 
and (jS]), after calculation with isospin variables one obtains the iViV-FSI amplitude in 
the form 

Mb = / [A^Mo + AgMi] (case 1), 

( 6 ) 

Mb = / [H^Mo — AgMi] (case 2), 

where /[••■] denotes the integration over the intermediate momentum (and sum over inter¬ 
mediate spin states) in the triangle diagram Mj, in Fig. [H Eqs. ([6]) show that the isovector 
and isoscalar yiV amplitudes Ay and Ag couple with 1=0 and 1=1 pn amplitudes Mq and 
Ml, respectively. This follows from the simple consideration that for the hnal TT^{pn) sys¬ 
tem with total 1= 0 (1) the pn system should be in an I = 1 (0) state. From Eq. ([S]), one 
can get “elastic” and “charge-exchange” pn amplitudes as 

M;' = M(pn^pn) = ^(Mi + Mq), M^" = M(np^pn) = ^(Mi - Mq), (7) 

and rewrite Eq. ([6]) in the more illustrative form 

Mb = / [{Ay+Ag)Mf^ - {Ay-Ag)M^l^] (case 1), 

, ( 8 ) 

Mb = / [{Ay-Ag)M^y^ - (A + H,)M“"] (case 2). 

The amplitude Mb is represented by two terms. In case 1, the hrst (second) term contains 
the photoproduction amplitude yp -X- vr^p (yn 7r°n) and subsequent rescattering of the 
fast proton (neutron) on the slow neutron (proton). Thus, the hrst term contains the 
elastic pn amplitude M®^ (fast proton —)■ fast proton), while the second - charge-exchange 
one Mp®* involves the fast neutron —>■ fast proton scattering. The changes for case 2 are 
obvious. In both cases, the relative sign “—” between two terms in Eq. ([8]) arises from 
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the isospin antisymmetry of the DWF with respect to the nucleons. In both Eqs. ([6]) 
and ([8]), we again follow only the isospin structure of the amplitude and show that the 
only difference between the results comes from the different relative sign of the term ^4^ in 
cases 1 and 2. Considering the vriV-FSI diagrams M^i and Mc 2 in Fig. [H we also obtain 

for the general case that both these amplitudes change from case 1 to case 2 due to the 
different sign of the terms, containing the isoscalar yiV —)■ ttN amplitude Ag. Hereafter 
in the present paper, we neglect the ttA^-FSI amplitudes, since their role was found to 
be small in Refs. Hi], and negligible at energies E above 200 MeV [5|. Then the total 
yd —)■ TT^pn amplitude can be expressed in the form 

M'yd = Mai + ^5 A = Ma2 + Mf), (9) 

where the lA amplitude Mai is the leading-order term (quasi-free production), and A 
includes the suppressed lA diagram Ma 2 and the AW-FSI diagram Mb. In quasi-free 
kinematics with fast and slow hnal-state nucleons, A is a relatively small correction term 
in comparison to the main one Mai- In this kinematics the relative momentum between the 
two hnal-state nucleons is large, and the contribution from AdV-FSI is suppressed, i.e., the 
Migdal-Watson effect [7| from s-wave AW-FSI at small relative momenta is suppressed. The 
quasi-free kinematics also implies not very small momentum transfer from the initial photon 
to the hnal pion, and thus, the FSI effect doesn’t essentially cancell the lA amplitudes due 
to the orthogonality of the hnal pn state to the wave function of the initial deuteron. We 
shall return to this point below in Sects. 4 and 5. 


3. Extraction Procedure for jN —>■ tz^N Reactions 

We now discuss the procedure of extracting the yA-reaction cross section from the 
deuteron data. First, consider case 1 (fast hnal-state proton, slow neutron). In quasi-free 
kinematics one has the relation [8] (more details are given in the Appendix) 


da 

dVt 


(yp Ti^p) 


1 da{jd) ^ 


E. 


p(l9.l)' 


( 10 ) 


Here: dVl is the solid angle element of the outgoing 7r° in the Tz^pj center-of-mass frame 
with the z-axis directed along the photon beam; E.y and are the photon laboratory 
energies in the yd process and the yp —)• 7r°p subprocess, and E'^/E^ = 1-|- {\qg\/Es) cosOg, 
where Eg and cos6g are the total energy and polar angle (z-axis along the photon beam) 
of the hnal slow neutron Ug in the laboratory system. The momentum distribution in the 
deuteron is given by p(g) = { 27 i)~^[u‘^{q) +w^{q)] (g = |g|) where u{q) and w{q) are the 
S- and A-wave parts of the DWF. It is normalized to Jp{q) dq = 1. One can improve 
Eq. fllOp . taking FSI corrections into account. With these corrections, we have 




Rp da {'yd) 


( 11 ) 


n{qg) dVLdqg 

The extraction procedure for the yp —)■ n^p reaction implies Eq. oa to be utilized with 
the deuteron cross section da{'yd)/dfldqg, taken from the data. Here, the correction factor 
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Figure 2: (Color online) The differential cross sections of the jp —)■ n^p (red solid curves) 
and 7 n —)■ (blue dashed curves) reactions at several photon energies = 340 (a), 630 

(b) and 787 MeV (c), which correspond to A(1232)3/2+, A(1440)l/2+, and iV(1535)l/2“ 
regions, respectively. 


Rp can be obtained from the model as 

_ da^P\jd) / dajjd) 

^ dfl dq^ / dVt dq^ ’ 

where the cross sections da^^\'jd) is calculated from the main lA diagram Mai, while da{'jd) 
is obtained from the full amplitude = Mai + A in Eqs. (jO]) and should restore the 
experimental data. Consider now case 2 (fast neutron, slow proton). Instead of Eqs. (fTTj) 
and flT^ . we obtain 


da Q Rn da{'jd) da^'^\'yd) /da{'jd) 

d^y'y^^- n{q^) dndq,' d^dq, / d^dq,' ^ ^ 

where dQ is the solid angle element of the outgoing in the vr°n/ center-of-mass frame 
with the z-axis directed along the photon beam. All other notations are the same as in 
Eqs. oa and (IT^ with the replacements p/ —)-n/ and Ug Ps- Here, the cross section 
da^'^\'yd) is also calculated with the main diagram Mai, but the fast nucleon is a neutron. 
The correction factors Rp in Eq. flT^ (case 1) and Rn in Eqs. flT^ (case 2) are calculated 
through the amplitudes given by Eqs. (0]) and ([6]). Going from case 1 to case 2, we change 
only the relative sign of the isoscalar photoproduction amplitude Ag in these Eqs. Thus, 
in the general case 

Rn ^ Rp- (14) 

However, when = 0 or = 0 in Eq. (|2]), we get = Rp. In this particular cases, the 
differential cross sections of the yp —?• 7r°p and yn reactions are equal. 

3.1 Comment on the A(1232)3/2+ Region 

Consider the 7r° photoproduction in the A(1232)3/2+ region. Supposing that in this 
region {E.y ^ 340 MeV) the yA —)■ vr^A proceeds through A(1232)3/2+, we obtain A^ = 0 
due to isospin conservation. In this case, we get Rn= Rp, i-e., the correction factors for 
the reactions yp —)■ 7r°p and yn —)■ vr^n appear to be the same, and we should also have 
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equal differential cross sections for these reactions. To illustrate this, we reconstruct in 
Fig. [2] these cross sections at several energies, obtained with the jN —)■ tiN amplitude 
taken from the SAID database [3]. Fig. [2]shows that proton and neutron cross sections are 
very close to each other in the A(1232)3/2’'' region {E^ = 340 MeV). At higher energies, 
the contributions from A(1440)l/2’'' and A(1535)l/2“ become important, the isoscalar 
amplitude Ag ^ 0 , and the difference between proton and neutron differential cross sections 
is visible. The main {Mai) and F SI-correct ion (A) terms in the yd —)■ n^pn amplitude 
M.yd dH]) also become different in cases 1 and 2 , and we get Ineq. flTTD . 

3.2 Comment on the ? 7 -Photoproduction Case 

Let us comment on the case of r] photoproduction on the nucleon {'jN rjN) and 
deuteron (yd —)■ ppn). The amplitude on the nucleon has isospin structure A{'yN —)■ 
rjN) = AyT 3 + Ag, which gives 

A(yp pp) = Ay + A^, A(yn pn) = -Ay + A^, (15) 

The difference with the vr^ case is that the isovector amplitude Ay changes sign for p 
production. Thus, the amplitudes Mai, Ma 2 , and Mf, on the deuteron with the vr^ replaced 
by an p are given by the same Eqs. (jl]), ([H]) and (E]) with small modihcations, i.e., the 
amplitude Ay (instead of A^) changes sign when going from case 1 (fast proton) to case 2 
(fast neutron). Thus, we arrive at the same conclusion (Eq. ITT)) in the general case 
of A^ 7 ^ 0 and A^ 7 ^ 0. Addition of pN-FSl amplitudes, similar to Md and Mc 2 in 
Fig. [H or the yd —)• ppn amplitude with intermediate pion photoproduction followed by 
the ttA —)■ pN subprocess also lead in general to Rp 7 ^ However, one should note 
that in this case FSI effects seem to be much smaller than for 7 r° photoproduction. The 
experimental data [141115] show that ~ 1 over a wide range of incident photon energies 
and p c.m. polar angles. A major difference to 7 r° photoproduction is that this reaction 
is dominated for not too high incident photon energies by the excitation of Fn resonances 
via the Eq+ spin-flip multipole. This means that in quasi-free production off the deuteron, 
the hnal state nucleons are dominantly in a spin-singlet, while for vr® production the spin- 
triplet configuration of the deuteron dominates the final state. This effect (different spin 
configuration) of initial (deuteron) and final np-system suppresses the role of FSI in the p 
case in comparison with the 7 r°-production case. The incoherent y-photoproduction off the 
deuteron was studied in a number of papers (see Ref. [12] and references therein), and the 
FSI effects were found to be important only in the very nearthreshold region. At higher 
energies, they are quite small [ 12 ], except for special kinematic configurations, where the 
lA contributions are suppressed. Theoretical results of Ref. |T2] show visible sensitivity of 
some polarization observables to pN and intermediate-pion {tiN —)■ pN) FSI effects, but 
this issue is beyond the scope of the present paper. 

4. Numerical Results 

Here, we perform some numerical comparison of the correction factors Rp fll 2 p and 
Rn UTS]) to illustrate Ineq. dn. Let us briefly describe the ingredients used in the calcu¬ 
lation of the reaction amplitude M^d dl])- We use the elementary pN^nN amplitudes. 
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generated through the GW pion photoproduction multipoles [10]. For the iViV-FSI term 
Mh in Fig. [H we include the s-wave pn-scattering amplitudes Mq and Mi, introduced in 
Eq. dSj). These invariant amplitudes for NiN 2 ^read 


Mo = 87rv5/^°)(p)((p+<Tp^)((p=+<T(pJ, Mi = 87ry/^f^l^{p){ip+ipl){ip^+ip^) (16) 


is the effective AW-system mass). Here: <^1234 are the nucleon Pauli spinors = 

1); = a 2 (p*. The amplitude Mo(Mi) with isospin / = 0(1) corresponds to the total spin 

S = 1(0). The amplitudes fpn^\p), where p = |p| is the relative momentum in the NN 
system, are expressed through the scattering lengths Oq 1 and effective radii Tqi, i.e.. 


f( 0 ,l) 
J pn 


{p) 


1 

-ag} + lrQ^iP‘^-ip' 



P^+^h 

q^+/3li 



(17) 


where {q,p) are the off-shell iViV-amplitudes, and q = \q\ is the relative momen¬ 

tum of the intermediate nucleons. The off-shell dependence of fpn^^°'^^{q,p) is introduced 
through the monopole-type formfactor, as given in Eq. fll7p . with parameters /9o,i- We take 
the typical value /9o,i = 1.2 fm“^ (the same as in Refs. [SIE|)- For oop and rop, we use the 
known values m- Oo = 5.4 fm, ro = 1.7 fm, oi = —24 fm, and ri = 2.7 fm. Hereafter in 
the present paper, we neglect higher partial waves (L> 0) in the iViV-scattering amplitude. 
To simplify the calculation of the iViV-FSI term Mi, [Fig. [1], we also take the amplitude of 
the jN —)■ ir^N subprocess out of the loop integral over the intermediate momentum in 
this term. This approximation allows us to calculate analytically the FSI term Mi,. The 
more precise calculations with higher NN partial waves included are planned for the future 
publication. The DWF was taken from the Bonn potential (full model) |13] . 

The results of the model for differential cross sections da/dfl* of the reaction yd n^pn 
for several photon laboratory energies E are compared in Fig. |3]to experimental data from 
Ref. [IT] [plots (a-/)] and to recent data from Ref. [18] [plots {g-m)]. Here, 6* is the polar 
angle of the outgoing in the c.m. system of the incident photon and a nucleon at rest with 
z-axis directed along the photon momentum. The dotted curves show the contributions 
from the lA term M^ = Mai-|- Ma 2 , while the dashed ones represent the results obtained 
with the full amplitude M^d= Mai-|- Ma 2 + Mi, (|9]). Comparison of these curves shows an 
important role of the A7V-FSI, which essentially decrease the cross sections at low energies 
in line with other papers [IHS]- At higher energies [plots(j-m)], the FSI effect is very small, 
except the region of small angles 6*, where the effect is sizeable. Our model predictions 
(dashed curves) sizeably overestimate the data, but the shape of the differential cross 
sections in the main is reproduced. It is quite possible that the full account of multiple¬ 
scattering diagrams (not only Mi, term in Fig. [1]) may lead to destructive interference and 
decrease the cross sections, insignihcantly affecting the shape of the distributions. 

Let us try to improve the theoretical description, taking into account that the 'yN —>■ ttN 
amplitudes in the diagrams M^i, Ma 2 and Mi, [Fig. [1] are not on the mass shell. Let us 
multiply the elementary jN —)■ irN amplitude by the off-shell correction factor of the form 


Fiq^r, q'^) 


A^ -F 

A 2 + g; 2 ’ 


(18) 


7 





where (g^) is the photon c.m. momentum in the 'fN —)■ ttN reaction, calculated with 
on-shell (off-shell) initial nucleon at a given effective 'yN mass. We add this factor to the 
lA amplitudes Mai and Ma 2 - For simplicity, we neglect this correction in the AW-FSI term, 
where the momenta of the nucleons in the deuteron vertex are effectively small in the loop 
integral, and the off-shell effect is expected to be small. On the other hand, we expect that 
the lA contributions at high energies [Fig,[3l plots (j-m), dotted curves], where FSI effects 
are very small, are overestimated in the model and should be suppressed. Solid curves in 
Fig. E] show the off-shell corrected results from the full amplitude dSD with A = 1 fm“^ 
in Eq. flTS]) (we use the value of the same order of magnitude as (3's in Eqs. ffT7|) ). The 
role of the off-shell correction is negligible at low energies [Fig. El plots (a-c)], but visibly 
decreases the cross sections at higher energies, and the theoretical description looks better. 
However, the model still overestimates the data, and varying the value of A fflSj) . one 
can not decrease the predicted cross sections in accordance with the data. Here we finish 
the comparison of the model with the data, leaving further developments, which perhaps 
should involve the multiple-scattering diagrams, to the next publications. 

Now consider the theoretical predictions for the correction factors Rp and i?„, intro¬ 
duced in Sect. 3. They are given in Fig. 0] at the same photon energies as in Fig. |2l For 
illustrative purposes the results are “averaged” over the spectator momentum q^, where 
Rp^n are dehned as 

^ _ da^P'^i'yd) /da{'yd) ^ _ da^^\-fd) /da{'yd) 

dn / dn ' dn / dn ' ^ ’ 

where the differential cross sections are integrated over the kinematic region \qj\ > |qr^|. 
Hereafter, we do not apply the off-shell correction Eq. flT8|) to the jN —)■ tiN amplitude. 
The plots at the top [(a-c)] show the factor Rp flT^ for case 1 with fast protons and 
slow neutrons, while the plots at bottom [{d-f)] show the factor for case 2 with fast 
neutrons and slow protons. The polar angle 6* of the outgoing pion is dehned in the tt^p 
(case 1) or vr^n (case 2) c.m. frame. The types of curves in Fig. 0] specify the calculation 
of the cross sections daipfd)/dVL, i.e., the denominators in Eqs. ([19]). The dashed curves 
are the results obtained with the full lA amplitude Mai + Ma 2 for these cross sections, 
i.e., they include the correction for the “suppressed” lA term Ma 2 - Addition of the NN- 
FSI term M;, leads to the solid curves. We see that the “suppressed” lA term Ma 2 alone 
already produces a sizeable deviation Rp^n ^ 1, which increases to small angles. The NN- 
FSI term M^, when included, considerably affects the results. Both terms {Ma 2 and Mb) 
essentially affect the results at small angles 6* (pions emitted at forward angles), where the 
conhguration with small relative momenta between the hnal-state nucleons dominates, and 
this 6* region narrows with the increasing photon energy. At higher angles, both effects 
are negligible, and Rp^n ~ 1- At = 340 MeV [plots {a,d)], we obtain a large effect 
at 9* ~ 0, where Rp^n ~ 6. We qualitatively interpret this results as a self-cancellation 
effect in the full amplitude M^d = Mai+ Ma 2 + Mb due to the non-ortogonality between the 
initial deuteron state and the hnal NN plane-wave state. This cancellation should enhance 
at small momentum transfers A from the initial photon to the final pion, decreasing the 
denominators da{'jd)/dQ in Eqs. flT^ . The effective value of A is minimal at = 0 






C 

"D 

\ 

b 

XI 


JD 

3 

* 

C. 

X 

b 

X 


4 

3 

2 

1 

0 

60 

50 

40 

30 

20 

10 

0 

12 




X 

b; 

X 

\ 

b 

X 



X 

3 

c 

X 

\ 

b 

X 



Figure 3: The differential cross sections of the 'jd n^pn for several values of the photon- 
beam laboratory energy E vs. 9* (the angle 9* is dehned in the text). The curves show 
the contributions: dotted - from the lA amplitudes Mai^a 2 [Fig. [I]; dashed - from the total 
amplitude dH]) with the AW-FSI. Solid curves include also the off-shell correction flTS]) 
for the 'yN -^tiN amplitude. The hlled circles: in the plots (a-/) - the data from Ref. [T7] 
(error bars include statistical uncertainties only); in the plots {g-m) - the data from Ref. [T8] 
(error bars include statistical and systematic uncertainties in quadratures), multiplied by 
1/A [A = 2 for the deuteron). 
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Figure 4: The correction factors Rp [plots (a-c)] and Rn [plots (d-/)], calculated according 
to Eqs. flT^ from the reactions 'jd ^ and ■yd ^ n^njPs, respectively. Left, middle 

and right plots - the results for = 340, 630, and 787 MeV, respectively. The numerators 
da^P’^^ydO. in Eq. flT^ are obtained from the leading lA amplitude Mai- Successive addition 
of the “suppressed” lA term Ma 2 and A7V-FSI term M^, when calculating the denominators 
da/dfl in Eqs. flT^ . leads to dashed and solid curves, respectively. 
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Figure 5: Comparison of the correction factors Rp (solid curves) and Rn (dashed curves) 
at the same energies E as in Fig. 01 Upper [(a-c)] and lower [(d-/)] plots correspond to the 
variants shown by dashed and solid curves in Fig. 01 respectively (see notations therein). 
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and decreases with the increasing photon energy. We shall also mark this point in the 
Conclusion. 

The same results are shown in Fig. E] as a comparison of Rp (solid) and (dashed) 
factors. The top [(a-c)] and and bottom [(d-/)] rows show these factors for successive 
inclusion of the “suppressed” lA term Ma 2 and NN-FSl term Mb for the denominators 
da{^d)/dVL in Eqs. ([19]). At = 340 MeV, i.e., in the A(1232)3/2+ region, we have 
Rp = Rn to a. good accuracy, since the role of the isoscalar 7 A —)■ nN amplitude Ag is 
negligible. At higher energies (630 and 787 MeV), we observe that Rpy^Rn, because both 
isovector and isoscalar amplitudes {A^ and A^) contribute. This difference is considerable 
at small angles 6* where the role of the “suppressed” lA term Ma 2 and A7V-FSI is enhanced. 
Thus, the results in Fig. [5] illustrate our analytical considerations discussed above, which 
lead to the general in Ineq. flTT|) . 


5. Conclusion 

We see that the values Rp^n strongly differ from 1 at small angles 6*. One obtains 
the large values Rp^n when including the A7V-FSI diagram, which decreases the differential 
cross sections at small angles. It means that at small momentum transfers the FSI term 
essentially cancels the lA terms due to orthogonality of the hnal pn-state to the deuteron 
wave function. This cancellation suppresses the denominators in Eqs. flT^ and (IT^ . i.e., 
the cross sections on the deuteron, and leads to the growth of Rp^n- We plan to return 
to this point with more details in the next papers. This effect, not connected with the 
Migdal-Watson effect [7] , was considered earlier in Ref. [IS] (see also Ref. II])- Concerning 
the physical meaning of the results one should remember that the photon has no dehnite 
isospin. For any given process, it may be represented as a superposition of I = 0 and 1 = 1 
states, i.e., I 7 ) = Co|0) + Ci|l), where the coefficients Coy depend on the reaction amplitude. 
Let us apply the charge symmetry (CS) operator [2D], which is the rotation by 180° about 
the y-axis in isospin space with the z-axis related to the charge, to the process 'yd^Tr^pfUs- 
We denote fast and slow protons (neutrons) as pf^s (^/,s)- The CS rotation of the initial 
and hnal particles gives 

I 7 ) ^ lY) = Co|0)-Ci|l) 7 ^ I 7 ), \d)^\d), |7r°) ^-|vr°), 

Ip)-)• |n), |n) ^ -|p), \pn)i^±\pn)i (+, - for /= 0,1). 

Then for the amplitudes of interest, we obtain 

M{-yd^7i^Pfns) = M{'y'd^Ti^nfPs) Y M{'yd^7i^nfPs). 

Here the hrst equality is the result of the CS rotation. The second inequality means 
that the corresponding amplitudes are different because the initial photons are different, 
i.e., I 7 ) 7 ^ I 7 '). As a result, we obtain Miyyd —>■ tt^p/Hs) Y Miyyd —>■ vr^ri/Ps), and that 
is the general reason for the Eq. ffTT|) . In the A(1232)3/2+ range, where the isoscalar 
7 A—)- 7 rA amplitude Ag is negligibly small, the proton and neutron correction factors are 
equal, Rp = Rn, to a good accuracy. This means that in this region one can extract the 
'yn^n^n differential cross section from the 'jd^ir^nfPs data by making use of the proton 


( 20 ) 

( 21 ) 
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factor Rp. This is convenient, since Rp can be directly obtained from Eq. CD. making 
use of the data on the free proton (yp —)■ vr^p) and on the deuteron (yd —>■ n^pfUg). In 
this way one can also verify the model predictions for Rp obtained from Eq. fll2p . In the 
region above the A(1232)3/2+, where both isovector and isoscalar yiV —)■ vriV amplitudes 
are important, we have Ineq. flTT|) . The main Rp/Rn differences in our results are observed 
at small angles d*, where contributions of the “suppressed” lA amplitude Ma 2 or AdV-FSI 
terms are important. However, there is a wide range of large angles 9*, where we have 
approximately Rp = Rn- To the end of this paper it is worth saying that the results of 
calculations, presented in Figs. EHSl were obtained not with the full program, but include 
a number of approximations with only s-wave AdV-rescatterings taken into account and 
simplihed version of the loop integral in the ESI term as was mentioned in Sect. 4. We 
leave more precise full calculations for the future publication. 
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Appendix: 


On the Relation between the Cross Sections on the Nucleon and Deuteron 

In quasi-free kinematics, the reaction jd^Tr^pfUg goes predominantly via the diagram 
Mai in Fig- E For definiteness, we consider case 1 with a fast proton and a slow nentron. 
In this approximation, the differential cross section on the denteron is related to that on 
the proton (yp vr^p) throngh Eq. ffTOj) . which is well-known [8]. Here, we give some 
explanations, which may be nsefnl when analysing experimental data. Hereafter, we use 
the notations 


= (E^, qr^), = (rUrf, 0), k^ = (E^, qr^), kj = {Ef, qr^), kg = {Eg, qrj (A.l) 

for the 4-momenta of the photon, deuteron, pion, fast proton and slow neutron (spectator), 
respectively. The total energies and 3-momenta in the laboratory frame (deuteron rest 
frame with z-axis along the photon beam) are given in brackets, and = \q^\ = E^, 
Qt, = \q^\ and = \qf^g\', 'rrid, nip, nin are the deuteron, pion, proton, and neutron 
masses. In the unpolarized case one can rewrite Eq. (nnii with the full set of variables as 


da. 


IP 


dfl 


{W, z-,m'p) = 


da. 


'yd 


E' 


n{qg) dVLdq^ 


{E^, z, p, qg, Zg), n(qrj = p(q,). 


(A.2) 


In Eq. (A.2) is dkl = dzdp {z = cos 6) the solid angle element of outgoing the 7r° in the 
7r°p center-of-mass frame with the z-axis along the photon beam, where 6 and ip are the 
corresponding polar and azimuthal angles. The laboratory polar angle of the neutron 
spectator is 9g with Zs = cos 9g. W is the effective vr^p mass related to qr^ by 


W‘^= {k,^+Pa-pgy = {E^ + TUd-Egf-E^-ql+2E,yqgZg {Eg = ^ml+ q^) (A.3) 

The momentum-distribution function p{q) and photon laboratory energy E!{ in the sub¬ 
process 7 p —)■ TT^p are defined in connection with Eq. jia. One can also write E!y = 

(IT^- mp)/2mp. 

For a given E,y, the cross section da,yd/dQ dq^ with unpolarized particles depends in 
general case on four variables, chosen in Eq. (A.2) as z, p, qg and Zg. The cross section 
da,yp/dQ, being a function of W and z, may also depend on the virtuality of the 
“target” proton in the subrocess yp —)■ 7r°p, where = {pd—PsY = {'^d — EsY — dl 7^ 
Thus, in view of Eq. (A.3) the cross section da,yp/dQ in Eq. (A.2) depends only on three 
variables, i.e., qz, Zg and z. Thus, the cross section da,yd/dfl dq^ on the deuteron should 
not depend on p in the model based on the leading lA diagram Mai in Fig. E 

For completeness, let ns express the cosine 2 ; in Eq. (A.2) through the particle momenta 
from Eq. (A.l). It can be obtained through the relations 


{kk)=E;{E:-q:z), 


^7 = 




-m: 


12 


2W 


Ef 


W^+m^ — nip 

2iy 


= (A.4) 


where E^ and E)) are the total energies of the photon and pion in the vr^p center-of-mass 
frame. 
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Let us rewrite Eq. (A.2) in a form more convenient for applications. Making the sub¬ 
stitutions dQ = dzdip and dqg = q^dqsdzgdips and averaging Eq. (A.2) over ip and psi we 
obtain 


da, 


IP 


dz 




d^a,yd 


27rn{qg) q^dqgdzdz. 


■ z, qz 5 Zs 


(A.5) 


In the real data analysis, one may divide the phase space of the reaction yd —)■ Ti^pn 
into small cubes, bounded by the values 


1a 1a 1a 

^ i 2 ^ 2 ^ 2 


(A.6) 


Let Aa{E^, z,qs, Zs) be the yd ir^pn cross section in such a cube (A.6). Then, we 


get 


da. 


7P 


dz 


{W,z;m'^^) = 


27rn{q, 


Aa{E^,z,qs,Zs 

AzAz^Aq^ q^ 


n{qs)= ^Pids), 

Jll/'y 


(A.7) 


where W'^ = W‘^{E.y,qs, Zg) =Eq. (A.3), E!y= {W‘^ — mp)/2mp, and m'p = {md —Eg)"^ — q^. 
At small qg (<^ mp^n), one may neglect the gf terms in m'p and (A.3). Then nip = nip, 
but W depends on and through the linear term 2E.yqsZs in Eq. (A.3), and E!) = 
E^[l -h {qg/Eg)zg]. 
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